Interaction of Lamb modes with two-level systems in amorphous nanoscopic 

membranes 
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Using a generalized model of interaction between a two-level system (TLS) and an arbitrary de- 
formation of the material, we calculate the interaction of Lamb modes with TLSs in amorphous 
nanoscopic membranes. We compare the mean free paths of the Lamb modes with different sym- 
metries and calculate the heat conductivity n. In the limit of an infinitely wide membrane, the 
heat conductivity is divergent. Nevertheless, the finite size of the membrane imposes a lower cut-off 
for the phonons frequencies, which leads to the temperature dependence ft oc T[a -\- 61nT). This 
temperature dependence is a hallmark of the TLS-limited heat conductance at low temperature. 

PACS numbers: 



I. INTRODUCTION 

The development of nanodetectors and the strict re- 
quirements on their performance triggered intense ex- 
perimental and theoretical studies of their thermal prop- 
erties. These detectors work usually in a temperature 
range around 1 K or below and are supported by thin, 
insulating membranes. The thickness of such membranes 
is of the order of 100 nm, which, in the given tempera- 
ture range, makes it comparable to the dominant ther- 
mal phonon wavelength. In problems where the phonon 
wavelength is comparable to or longer than some of the 
dimensions of the system in question, the three dimen- 
sional (3D) phonon gas model cannot be applied anymore 
to calculate the system's thermal properties. Instead, 
one has to use the phonon modes specific to the system, 
which are the eigenmodes of the elastic equation for the 
given geometry. 

The membranes that support the detectors are made 
of amorphous, low stress silicon-nitride (SiNx) and their 
thermal properties have been measured in various ge- 
ometries by different groups (see for example Refs. 
[lii2n3i,4,i5n6|] ) . Depending on the quality and the dimen- 
sions of the samples, and possibly the temperature range 
in which measurements were done, the heat flux along the 
membrane may be due to either diffusive^i^i^ii or radia- 
tive phonon transport In the case of diffusive phonon 
transport, it was observed that the heat conductivity k 
is roughly proportional to . 

For a better thermal insulation of the detector, the 
underlying membrane is sometimes cut. The result is a 
self-supporting structure with a wider area in the middle, 
connected to the bulk by long, narrow bridges, like in 
Fig. [T]3fii^ The heat conductivity along such bridges has 



FIG. 1: (a) Full, dielectric membrane, (b) Cut membrane for 
better thermal insulation. 



again a power law dependence on the temperature, k cx 
TP, where p takes values between 1.5 and 2.— 1^ For the 
samples measured in p increased with the width of 
the bridge. 

The heat capacity of a membrane is more difficult to 
measure directly, since the membrane is always in contact 
with the bulk. However, it can be estimated by applying 
AC heating and measuring the amplitude of the tempera- 
ture oscillations. In this way Leivo and Pekolai^ observed 
that the ratio cy/K, where cv is the heat capacity, in- 
creases with temperature for the narrowest bridges. 

As mentioned above, to explain theoretically all these 
observations we have to work with the proper set of 
phonon modes. For wide membranes with parallel sur- 
faces, the eigenmodes of the elastic equations are called 
Lamb modes and horizontal shear modes, as explained 
for example in Ref. Q . Using these modes and their dis- 
persion relations, we could describe quite well the ther- 
mal properties of the membranes in the low temperature 
limit^ i.e. at temperatures where the the characteristic 
thermal wavelength of the phonons is much longer than 
the thickness of the membrane. 

Nevertheless, the same temperature dependence of 
the heat conductivity and heat capacity persists also 
in a temperature range where this characteristic thcr- 
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mal wavelength becomes equal, or even smaller than the 
thickness. This can no longer be explained only by the 
dispersion relation of the Lamb modes in the membrane 
and we have to take into account the amorphous struc- 
ture of the membrane and the resulting phonon scatter- 
ing. 

An amorphous material contains dynamic defects 
which can be modeled by an ensemble of two-level sys- 
tems (TLS)ji2iiiiii^ A TLS can be understood as an atom 
or group of atoms which can tunnel between two close 
minima in the configuration space. Any deformation of 
the material disturbs the TLS, which can have a transi- 
tion (an excitation or a de-excitation) to the other en- 
ergy level. A passing phonon produces such a deforma- 
tion and therefore may be scattered by the TLS. In bulk 
materials, the phonon modes are simple, trans versally 
or longitudinally polarized plane waves and the defor- 
mation field they produce can be described by only two 
parameters, the wave vector (or the wavelength) and the 
polarization. As a consequence, in the so called standard 
tunneling model, the interaction Hamiltonian has a very 
simple structure . In a mesoscopic system, the de- 
formation caused by the displacement field of a phonon 
mode is more complex and the TLS-phonon interaction 
Hamiltonian has to be modified accordingly. This was 
done in Ref. [isj . Here we use this more general Hamil- 
tonian to calculate the interaction of the phonon modes 
of the membrane with the TLSs. 



II. TLS-PHONON INTERACTION 

A TLS is described by a Hamiltonian which has the 
form 



H' 



TLS 



A 



A 



(1) 



when written in the basis formed by the ground states 
of the two potential wells between which the system 
tunnelsi^. In Eq. (IT|) A is the asymmetry of the potential, 
A is the tunnel splitting, and ax, (Jz are Pauli matrices. 
The Hamiltonian ^ can be diagonalized by an orthogo- 
nal transformation O, 



H' 



TLS 



(2) 



and we obtain the excitation energy, e = VA^ + A^. Ev- 
erywhere in this paper the superscript T denotes the 
transpose of a matrix. 

The TLS parameters A and A are not the same for all 
the TLSs in the material, but they are well modeled by 
the distribution P(A, A) = Pq/A in the unit volume of 
the material. We can rewrite the function P in terms of 
the more practical variables e and u = A/e, 



Pie,u) 



Po 



(3) 



The TLSs that have an excitation energy comparable to 
fceT are very efficient phonon scatterers. 

The deformation due to the displacement field of a 
phonon is quantitatively described by the strain field, 
which will be represented here by the 6-componcnt vector 
S.— If we denote the displacement field by u(r), then the 
strain is defined as the symmetric gradient of u(r), i.e. 

= (Vsu)^ = [dxUx,dyUy,dzUz,dyUz +dzUy,dxUz + 
dzUx, dxUy -\- dyUx)- This deformation adds a time- 
dependent perturbation to the Hamiltonian ([!]), which 
we shall denote by Hi. The perturbation is assumed to 
be diagonal, when written in the basis of the two ground 
states of the potential wells (hke in Eq. ^) fi2^1LlLl±l^ 



Hi 



(4) 



(6) 



and linear in the strain field at the location of the 
TLS^^ 

S = 27T* • [r] • S . (5) 

The other quantities in Eq. ([5]) are the coupling con- 
stant 7. the six component vector T, which is de- 
fined by a generic orientation t of the TLS as T = 
{tl,ty,tl,2tytz,2tztx,2txty)'^ , and the 6x6 matrix of 
the deformation potential parameters [r]. For isotropic 
materials, the matrix [r] is 

/ 1 C C \ 
C 1 C 

C C 1 

^00 

e 
Vooooo^y 

with the TLS potential parameters ^ and C that satisfy 
the condition ^ -|- 2^ = 1 

To calculate the scattering probabilities we have to 
write Hi in the second quantization. For this we de- 
note the TLS excited state by |t) and its ground state by 
I J,) and we introduce the "creation" and "annihilation" 
operators a^ and a, respectively, so that a^||) = |t), 
a|t) = a^lT) ~ 0; and a|J,) = 0. The operators a and 

obey Fermi commutation relations and in matrix form 
we have CTz = (2a^a — 1) and ax = {a^ + a). The bosonic 
creation and annihilation operators for phonons will be 
denoted by bj^ and 6^, respectively, where n stands in 
general for the quantum numbers of the phonon modes 
(see for example Rcfs. [lJI3)- Using these notations and 
applying the transformation O to the total Hamiltonian 
{O^IHtls + Hi)0 = iJ^LS +H{),wc obtain 



-^T^ . [r] ■ J2 [%b, + S;bl] (at + a) . (7) 



E 



[SA + S;6t] (2ata-l) 



In the first order perturbation theory, the phonon ab- 
sorption and emission rates are determined by the off- 
diagonal elements of H'l, i.e. of the second row of Eq. 
([7]). Higher order processes are not considered here. 
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A. The phonon modes in the membrane 

The phonon modes of a free standing, infinite mem- 
brane are divided into three groups, according to their 
symmetry properties. One group is formed of simple, 
transversally polarized modes, called the horizontal shear 
modes (h). The two other groups are the symmetric (s) 
and antisymmetric (a) Lamb modes. Together, these 
modes form a complete, orthonormal set of functions 
for the elastic displacement fields in the membrane and 
their proper quantization has been carried out in Ref. 
[l6l | . In this paper we shall use the results and notations 
from there and we shall call the three different types of 
phonons (i.e. h, s and a) polarizations. 

We assume that the membrane of thickness d is placed 
parallel to the (xy) plane and its parallel surfaces cut the 
z axis at ±d/2. The phonons propagate in the {xy) plane 
with the wave vector ky = fc||k|[, of real fcy. We use hat 
to denote unit vectors. 

Along the z direction, the phonon modes are station- 
ary. As the h modes are pure transversal waves, they 
have only one wave vector component along the z direc- 
tion, which we denote by kh- The s and a waves are 
superpositions of transversal and longitudinal waves of 
wave vector components along the z direction denoted 
by kt and ki, respectively. Due to the boundary con- 
ditions, which demand that the membrane surfaces are 
stress- free, kh takes the discrete values mn/d, with m 
taking all integer values between and oo, whereas kt 
and ki satisfy the more complicated relations^ 

tan(fctd/2) ^ 4fctfc;fc|| 

tan(fc,d/2) " - fc2)2 ' i^aj 

for the symmetric modes and 

tan(fc;d/2) ^ 4fctfc;fcj| 

tan(fctd/2) (fc2 - fc2)2 ■ ^ ) 



Uh 

Us 



Ua 



As one can see, when kt or ki take imaginary values, the 
trigonometric functions in ^ will switch into hyperbolic 



for the antisymmetric modes. Equations ([8]) plus Snell's 
law, ^2 = Cj (fcjj + kt) = c2(fc2 + kf), enable us to write 
ki as a function of kt for each of the polarizations s and 
a. (In Snell's law, q and q are the respective transversal 
and longitudinal sound velocities for bulk media.) Like in 
the case of the h modes, the dispersion relations for the 
symmetric and antisymmetric modes split into branches, 
i.e. Eqs. ([5]) and Snell's law do not give only one function 
kiikt) for either set of modes, but produce an infinite, 
countable set of such functions. These functions will be 
called phonon branches and we shall number them with 
771 = 0, 1, . . ., as we did in Ref. [l^, where branches of 
bigger m lie above branches of smaller m. 

A simple way to express the quantum numbers of 
the phonon modes is to use Eqs. ([8]) and Snell's law to 
write the hmctions fct(fc|[) and /s/(fcj[). Then each branch, 
of polarization a = s^a and branch number to, is go- 
ing to be described by the continuous set of numbers 
[fci(A:||), A:t(fc||)]o-,m, with /cy taking values from to oo. 
We therefore choose the set /i of quantum numbers that 
specify the phonon modes in Eq. ([71) to be /i = {cr, to, ky }. 

The functions fct(fc||) and ki{k\\) may take both, real 
and imaginary values. To distinguish between these sit- 
uations, we write the imaginary values of kt as int and 
the imaginary values of ki as ini. In these notations, ki^t 
and Ki^t take always positive, real values. 

In order to simplify the later discussion, we shall re- 
place kt and ki with the complex quantities kt = kt + int 
and ki = ki + ini, respectively. Note, however, that kt 
and ki are never really complex, but they are either real 
or imaginary, as long as /cy is realJ^ 

The displacement fields of the phonon modes are^il^ 



(9a) 



(9b) 



(9c) 

I 

functions. The normalization constants Nh, Ns and Na 
arc given hy^ 



I 

= Nh cos {kh{z - d/2)) (icy X i)e^(^r'r-^t) 

= Ns jifct (2fcy cos{ktd/2) cos(fc;z) -I- [kf - k^] cos(fc/d/2) cos(fct2:)^ ky 

- fc|| (2ktki cosiktd/2) sin(fc;z) - [kf - fcjj] cos(fc,d/2) sin(A:tz)) z| e*^''" 

= Na ^ikt (2k1 sm{ktd/2) sm{kiz) + [kf - kf^] sin(fc,d/2) sin(fct2)^ ky 

+ fcy (2ktki sin(fctd/2) cos(fc,z) - [kf - kf,] sin(fc,d/2) cos(A:tz)) z| g'^'^ii . 



4 



^-^'-{1,2 III (10^) 
-4fc{||cos(fc/d/2)p ^Ktdfctp + sinh(Ktd) - A:t(|fctp - /fc^) sin(fctd)) | (10b) 

-Akl\&in{hdl2)\'' i^nt{\h? + kl)s\nh{ntd) + kt^^^^^ \ , (10c) 



where A is the area of the membrane and V = d- Ais the 
volume. To obtain the expressions for Ng and Na for the 
different combinations of real or imaginary kt and fc; , one 
has to take the limit to of their redundant components 
in (jlOb|) and (|10cp . For example if kt is real and ki is 
imaginary, we calculate the corresponding normalization 
factor by taking in Eq. (jlObp or (jlOcp the limits 
and ki 0^ 

III. TRANSITION RATES 

Now we have all the ingredients to calculate TLS tran- 
sition rates or phonon absorption and emission probabil- 



I 

ities. We shall denote the phonon-TLS quantum states 
by \nfj,,l) or jn^,!), where we denoted the number of 
phonons on the mode /i by n^. Using Eq. ([7]) we write 
the emission amplitude of a phonon by a TLS as 



{n^A\Hi\n^ + lA) = -Vy2^^^'" ^^^^ 
with Mf^ given by 

Af^(t)=T^.[r].S^. 
Explicitly, for the three phonon polarizations we have 



Mh,„^,k^^{i) = 2S,Nh {-tyt,khs\n{kh{z - d/2)) + itJyk\\cos{kh{z + d/2))] e^'''^^'' (12a) 
Ms,n^,k,{i) = Ns [-2ktk\\ cos(fctd/2)cos(fczz){fcjj[C + (1 - C)i'] + kf[C + (1 - C)i']} 
+fctfcj| [kt - kj] cos{kid/2) cos(fctz)(C - l)(t^ - t^) 

-8it^t,^ktkikfcos{ktd/2)sm{kiz) - 2it^t,£,[k'f - cos(A:id/2) sin(fctz)| (12b) 
Ma,m,k, (t) = Na {-2fctfc|| sin(^td/2) sin(fc,z){fc|[<2 + ({I - tl)] + kf[tl + C(l - tl)]} 
-fc|,fct[fc2 _ sin(^zd/2) sin(fctz)(l - C){4 - tl) 

+8it^t^^ktkikj sin(fctd/2) cos(fc,z) + 2tJ;,^[kt - fc||]^ sin(A:/d/2) cos(fctz)| (12c) 

I 



where kt and ki are implicitly determined by the branch ^enn respectively, 
number, m, and fcy . Using Fermi's golden rule, we calcu- 
late the phonon absorption and emission rates T'^^,^ and -2*2 

-A' _ ^ ^ ^ I A T |2 



abs '^'^'^^ ~2 A 2 
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Tom = —^\M^Wn,. + l)5{hu:,. - e) , (13b) 

where e is the energy of the TLS, as defined in Section 
im and to is the angular frequency of the phonon. 

In an amorphous sohd, the orientations of the TLSs 
are arbitrary, so the relevant quantities for our calcula- 
tions are the averages of r^^j^^ over the directions t of the 
TLSs. The only quantity that depends on t in the Eqs. 
(fT3|) is Additionally, we assume that the distribu- 



tion of TLSs in the membrane is uniform, which leaves 
again |M^p the only quantity dependent on z in the ex- 
pressions for the absorbtion and emission rates. As we 
are interested in an average scattering probability rather 
than in a detailed description of where along the z di- 
rection the scattering takes place, we also average |M^p 
along z. Denoting by (•) the average over the TLS orien- 
tations and the z variable, we obtain 



{\Mh\') 

(|M.|'> = ^-^\^Ci\h\'kl\kl + kt\'\cos{ktdl2)\' 



V 



sinh(K;(i) sin(fc/(i) 



2ki 



2ki 



+Ct\~k'^~kff\cosCkid/2)\^ [{\kt 



-2Ctfct(2fcf-fc||fct|2- |fci|6)|cos(fcid/2)|2sin(fctd) 
-2CtKt(2fcf + fcf|fct|2 + |fct|6)|cos(fcid/2)psinh(Ktd) 



2h 



Af2 _ 

{\Ma\') = -2{^Ci\h\'kj\kj + kf\^\smCktd/2)\^ 



/ sinh(Kid) sin(fci(i) 



2ki 

2 2\Sinh(Kt(i) 



+Ct\kt - fc^n sin(fc,d/2)|^ + fc^) 

-2Ctfct(-2fcf + fc||fct|2 + |fc,|6)|sin(fc,d/2)|2sin(fctd) 
-2CtKti2kl + kflkif + |fct|^)|sin(fc,d/2)psinh(Ktd) 

I 



2k, 



ilhl^^klfJ^ 



(14a) 



(14b) 



(14c) 



Ct and Ci are constants that depend on the deformation 
potential parameter ^, Ct = 4^^/15 and Ci = (15 — 40^ -I- 
32^2)/i5^ Note that Q > C* > for all Using 
Eqs. (fT3|) and (fT4|) . we can calculate the TLS and phonon 
relaxation times. 



coth(/3e/2) J2 —{\M,,\^)Sihu;^^e) (15) 



p 



and 



~ 9 

n 7^ 



- — y u'(|A/mI'> tB.nh{f3e/2)5ihij^ - e) 
^^fXI^{\M^\^)tanh{Phuj/2), 



(16) 



respectively. Here we changed the summation over u and 
e into a two-dimensional integral and used the TLS den- 
sity ©. 



Using Eq. we can calculate the two-dimensional 
heat conductivity along the membrane. 



A 



dT 



(17) 



1 



( duj^y^ k\\ujl coth{l3huj/2) 



\dk\\) (|M^|2) sinh2(/3?icj/2) 



At an arbitrary temperature, n has to be calculated nu- 
merically. In this paper we give only the analytical low 
temperature approximation. 



IV. LOW ENERGY EXPANSION: 
ASYMPTOTIC RESULTS 

We can analytically calculate the scattering times or 
the thermal properties of the membrane only in the long 
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wavelength limit, i.e. for the branch m = for each of the 
three polarizations of the phonon modes and fcy <g; l/d. 
The calculation of thermal properties in this limit corre- 
sponds to a temperature range in which k^T <^ hct/d. 

First we have to calculate the relaxation times for each 
polarization. 



A. h mode 

For the lowest branch of the h mode, g is linear 
in fc|| and using Eq. (|14ap . the calculation of Th^o.k^^ is 
straightforward for any fcii. 



hpct 1 coth(/3?itj/2) 
tt^^Pq Ct TwJ 



(18) 




FIG. 2: The ratio rh^o.fey /Ta,o.k^^ as function of the ratio Ct /cf 
for different values of Ci/Ct- The value Ci/Ct ~ 4/3 is the 
minimum possible value and thus all other curves lie above 
the solid curve in the plot. 



s mode 



To get a long-wavelength expression of the dispersion 
relation for the lowest branch of the s mode, we note 
that ki ~ iKi takes imaginary values, which turns Eq. 
iil into 



tan(fctV2) 



tanh(KiV2) (k? - kfA 



AktKikj 
2T2 



(19) 



We expand the trigonometric functions in Eq. (jl9p to 
leading order and obtain 



<^s,0,fe|i 



(20) 



Using this, we calculate the relaxation time for this 
branch, 



T's,0,fc|i 



4c?(c? 



c?) coth(/3?iw/2) 



_ hpct 1 coth(/3?iw/2) 



2c?) 



(21) 



Nevertheless, this asymptotic expression is not enough 
for the calculation of (|Ma^o,fc|| P), as it turns out that 
both expression (|10cp and expression (|14cp are zero in 
the first and the second leading orders. Therefore, we 
have to expand to the third leading order to get non-zero 
results. From Eq. (|8bp we obtain 



2m* 



9Qcf 



(24) 



from which we finally get 



Ta,0,fc|i 



hpc^ 



cot'h{f3huj/2) 



nj^PoCict + Ctcficf-2ci) 
hpc^ 1 coth{phu;/2) 



huj 



TTj'^Po Ca hoj 
with Ca = 4Cs. 



(25) 



D. Comparison of the scattering rates and mean 
free paths 



a mode 



The antisymmetric modes have a more complicated 
asymptotic expansion. First let us remark that for the 
lowest branch and any fcy , both ki = ini and kf = ik^ 
take imaginary values, so we write Eq. (|8b|) in the form 



tanh(Ki(i/2) 



AKtKlku 



tanh(Ktc?/2) (kj 



(22) 



Expanding this equation to the second leading order, we 
obtain a quadratic dispersion relation for very small fcii ^ 



Wa,0,fc|| 




(23) 



The first thing to observe is that, although the disper- 
sion relations for the s and a modes are different in the 
low fcy limit, Ts,o,fc|| and o,fc|| are related by the simple 
equation Ts,o,fc|| = 4Ta,o,fc,| ■ In other words, the scatter- 
ing rate for the s phonons is 4 times smaller than the 
scattering rate of a phonons at the same to. 

Let us now compare Ta,o,fe|| with Thfl^kn ■ ^'^^ t^^-S, we 
calculate the ratio 



T;i,o,fc|| 



Qct + acncf^2c^,) 



A~[CilCt){ctlcif 
1 - (q/q)2 



(26) 



In any normal material (i.e. with positive Poisson ratio). 
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the ratio Cj/c^ is restricted to < c^/c^ < 1/2 and for 
Ci/Ct we haveH Ci/Ct = 15/4^2 „ ^q/^ + 8 > 4/3. In 
Fig.O we plot the ratio Th^.k^jTafl^k^^ as function of Cj /c^ 
for different values of Ci/Ct- 

We first remark that in the limit q/q 0, 
Tft,_o.fe|| /Ta^o,fc|| = I7 independent of the value for 
Ci/Ct- Increasing ct/ci will result in a decrease of 
T/i.o.feii /Ta.o.fcii , until a minimum is reached at {ct/ciY = 
1 - yjl - {Ci/Ct)-\ which hes between and 1/2. Af- 
terwards T/i^o,fc|| /■'"a,o.fe|| increases monotonically until it 
reaches the value {Ci/Ct)/2 for {ct/ciY = 1/2. As 
T^i.o.feii /Ta,o,fe|| increases monotonically with CijCt^ we 
conclude that Th.Q.ki^ < Ta,o,k^^ for any q/q, as long as 
Ci/Ct < 2. For Ci/Ct > 2, t/i^o.^h '"^'^ either smaller 
or greater than ra,o,fe|| i depending on whether {ct/ciY 
is smaller or greater than Ct/Ci, respectively. A typi- 
cal value for [ct/ciY in SiN^, is 0.36, which means that 
'Th,o,k\^ < ''a,o,fc|| as long as Ci/Ct is smaller than 2.78. 

When comparing Ts^o.fcn ^^^^^ o./cy , we encounter 
a similar situation. As Ts^o./ty = 4ra,o,fe|| the ra- 
tio T/i,o,fc|| /Ts,o,fc|| has the same features as the ratio 
''"/i.o.feii /''a,o,fc|| J except for the fact that the critical value 
for Ci/Ct is 8, i.e. T/j^o./cy is always smaller than Ts^o,fcy i 
if Ci/Ct < 8, and can be either smaller or greater 
than Ts_o.fe|| for Ci/Ct > 8. For the SiN^ typical value, 
[ct/ciY = 0.36, we have Th.o.kn < ^s.o.fcy as long as 
Ci/Ct < 17.6. 

More interesting than the scattering rates, is to com- 
pare the phonon mean free paths, since these can be di- 
rectly measured experimentally. For this, let us first use 
the dispersion relations (PO)) . (^5]) . and i^hfi^k^ = Qfeji to 
write the expressions for the mean free paths: 

fipcl 1 coth(/3ntj/2) 



■k^'^Pq Ct hw 
fipc] 1 coth(/3nctfc|t/2) 
7772 Po 'Ct fik\\ 



CsTs.O.ku 



hp4 2^/1 - cl/cf coih{(5huj/2) 
hpcj 1 coth(/3?iCsfc|[/2) 



(27a) 
(27b) 

(28a) 
(28b) 



la,0,l 



2nuj 



-Ta,0,k„ 



hpcl 2(l-c^/c^)V 4 
coth(/??itj/2) 



(29a) 



TTJ^Pq Ca hk\\ 



A way to determine the mean free path of phonons is to 
measure the resonant attenuation of ultra-sound, propa- 
gating along the membrane. If this is experimentally im- 
possible, another way to determine the material param- 
eters is to make acoustic measurements on thicker and 
wider membranes. Note however that elastic waves atten- 
uate not only because of resonant scattering of phonons 
(Eqs. HTJini), but also due to energy relaxatiouiii Never- 
theless, since we are interested here in the thermal prop- 
erties of the membranes, only the resonant scattering is 
important and we disregard the energy relaxation mech- 
anism. 

To analyze the results ^7^, ([28l), and ((291) in more de- 
tail, we expressed the mean free paths both in terms of 
the angular frequency and in terms of /cy . If the elas- 
tic modes of different polarizations are produced with 
the same to, then we should compare the mean free 
paths as given by the expressions (|27ap . (|28ap . and (|29ap . 
which we denote as lcr,o,ui- For example, lh,a,uj/ls.o,ui = 
iTh,o,uj/Ts,o,uj)-{ct/cs) which is smaller than (rhfi^uj /ts, o^uj) 
in any material. The discussion we made above about 
Th.o,uj/Ts,o,uj applies here too. 

The expressions ((29|) for Za.o.w are very different from 
the ones for Ihfi.cj and Isfi^Lj- Nevertheless, the expres- 
sions (1271), (12S1), and (1211) are calculated for k\\ < 1/d, 
which means that ud <C (27r/fcj|)cij = q which implies 

\fujdfct 1. Taking this into account when we compare 
the expressions for lh,o,uj, 's,o,w, and la,o,ujj at the same 
w, we conclude that, as a function of frequency, for low 
enough frequencies the antisymmetric Lamb modes have 
the shortest mean free path. This is a consequence of the 
fact that the group velocity of the a modes decreases to 
zero as fell decreases. 



If we compare I 



and I 



s,0,A;|| 



as functions of 



{29b) 



fcy, we see that /h,o,fcy/^s,o,fcy = (•^/i.c/ty /Ts,o,fe|| ) • 
[coth(/3?ictA:|[/2)/ coth(/37iCsA:y/2)], which is bigger than 
•r/i,o,fe||/Ts,o,fe||, since Ct < Cs implies coth(/3?ictfcy/2) > 
coth(/3?iCsfcj| /2). Comparing la,o,k^i with the expres- 
sions p7b[) and (|28bp . we observe, for example, that 
both Za,o,fe||/^/i,o,fe|| and /a,o,fe|| /^s,o,fey are proportional 
to coth(/3?itjQ^o,fc||/2)/coth(/37iij;ft,o,fcy/2). But again, for 
long wavelengths, due to the quadratic dependence of 
^a,o,fc|| on A:y, we have uja,o,k^i < '^h/s,o,k^^- Moreover, 
coth(a:;) ~ l/x for x 0, so both la,o,k^^/lh,o,k^^ and 
^a,o,fc|| /^s,o,fc|| are proportional to l/Zcy and become very 
big in the limit of long wavelengths. In conclusion, as 
function of fc|| in the limit dk^^ ^ 1, the antisymmetric 
modes have a much longer mean free path than the sym- 
metric and the horizontal shear modes with the same kn . 



E. Calculation of the heat conductivity 

Now we calculate the heat conductivity in the limit of 
low temperature. In that limit, only the lowest branch of 
each polarization will be occupied and we can write Eq. 
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(fT7|) in the form 



CONCLUSIONS 



167r kBT^ 



-y 



dio- 



smh^{(3huj/2) 



, (30) 



where the fowcr hmits uj* q were introduced for the rea- 
sons that will become clear imediately. Using Eqs. (|27t - 
I29p for the mean free paths we express k as a sum of 
three contributions: 



■T 



iKo) /(<o) 2/(x*,o) 



(31) 



where x*q = (3hLj* Q and by I{x) we denoted the integral 



lix) 



dy 



coth(ty/2) 
sinh^(2//2) 



2^x 



(e^ - 1)2 



1 



81n(l 



(32) 



Note that, although the mean free paths for the h and s 
modes have different functional dependences on u> than 
the mean free path for the a modes, the integrand in Eq. 
(PTjl is the same for all three modes. The role of the lower 
cut-off in Eq. ((30|) becomes obvious when we look at Eq. 
([32ll : the integral I{x) has a logarithmic divergence in 
x = 0. 

If the cut-off is small enough, then we can approximate 
I{x) by 



I{x) w 12- 81n(x) 
and inserting this into (|3ip we obtain 



(33) 



T 



3 - 2\n{f3nuj* „) 



3-21n(/3?i<o) , 2[3-21n(/3;i<o)] 



Ct 



Cs 



Ca 



,(34) 



where the first, the second, and the third terms in the 
square brackets above give the contributions of the h, 
s, and a phonon modes to the heat conductivity. The 
above expression leads to the temperature dependence 
K oc T(a -|- 61nT) and this dependence is a hallmark of 
the TLS-limited heat conductance at low temperature. 

For a numerical estimate let us use for the cut-off the 
finite size of the membrane, which limits the wave vec- 
tors to values of the order of ^TxjyfA. For the typical 
experimental parameters T = 0.1 K, \fA = 400 ^m 
and d = 200 nm^ii^i we have ln(x;i g 27r/L) ~ —4.9, 
ln(a;s,o,27r/L) ~ -4.4, and ln(a:a,o,27r/L) ~ -11.4. But 
since Ca = ACs (see Eqs. (PT|) and (^5]) ). the contribu- 
tions of all the phonon polarizations to the heat conduc- 
tivity arc of the same order. 



We used the model introduced in Ref. [lj| to calculate 
the scattering of the elastic modes in a thin, amorphous 
membrane. We modeled the scattering centers in the 
membrane by an ensemble of TLSs with the same prop- 
erties and distribution over energy splitting and asymme- 
try as the TLSs in a bulk material. If this assumption is 
valid remains to be checked by experiment. We obtained 
the expressions for the TLS relaxation time (|15p . for the 
phonon scattering time (|16p and for the heat conductivity 
K (dll). 

For general temperatures, the heat conductivity and 
the scattering times have to be calculated numerically. 
We calculated analytical low temperature approxima- 
tions and compared the mean free paths of different 
phonon polarizations. In this way we observed that the 
contribution of the lowest branches of the phonon modes 
to the heat conductivity are logarithmically divergent at 
/cjl 0. This could be a reason for which in some ex- 
periments a radiative heat transport is observedi^ Nev- 
ertheless, there is a natural lower cut-off of k\\ due 
to the finite size of the membrane. This cut-off renders 
K finite, which, in the low temperature limit behaves like 
K oc T(a + 61nT). This behavior is a hallmark of the 
TLS-limited heat conductance at low temperature. 

Due to the dispersion relations of the phonon modes, 
the TLSs distribution in the low energy limit has a big- 
ger impact on the heat conductivity in thin membranes 
than in bulk materials. If we for instance modify the 
distribution Q into 



P'{e,u) 



Pn 



with an extra energy dependence, e~", we make the ex- 
pression (j30p for K convergent even in the uj* q ^ limit, 
which leads to a low temperature asymptotic dependence 
of K oc But if this is the situation or not has to be 

decided experimentally. 
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